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ABSTRACT 

A  waterflood  of  the  Second  Wall  Creek  Sand  of  Naval  Petroleum 
Reserve  No.    3  is  being  considered  in  order  to  increase  the  ultimate  oil 
production.     A  streamline  model  of  the  Second  Wall  Creek  Sand  was 
developed  here  and  used  to  predict  the  response  of  the  reservoir  to  a 
waterflood. 

Results  of  the  model  study  indicate  that  the  oil  production  of  the 
Second  Wall  Creek  Sand  can  be  increased  by  as  much  as   15,  000,  000 
barrels  over  a  twelve  year  period  by  waterflooding.     This  would  bring 
the  total  production  available  from  the  Second  Wall  Creek  Sand  to  around 
32,  000,  000  barrels. 
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CHAPTER  I 
INTRODUCTION 

Naval  Petroleum  Reserve  No.    3  is  located  in  Natrona  County, 
about  thirty-five  miles  north  of  Casper,   Wyoming.     The  major  geo- 
logical structure  on  the  reserve  is  an  anticline  called  "Teapot  Dome". 
The  Second  Wall  Creek  Sand  is  the  most  important  of  the  several  oil 
bearing  formations  on  the  anticline. 

The  Second  Wall  Creek  Sand  is  a  solution  gas  drive  reservoir. 
There  appears  to  be  very  little  natural  water  drive.     Early  estimates, 
based  on  recovering  one -seventh  of  the  original  oil  in  place,    predicted 
the  recoverable  reserves  from  this  formation  to  be  around  17,000,000 
barrels(l). 

Waterfloods  of  the  same  formation  in  neighboring  fields  have 
proven  to  be  successful  in  increasing  production.     Therefore,    in  order 
to  enhance  the  recovery,   the  development  and  production  of  the  Second 
Wall  Creek  Sand  as  a  waterflood  is  being  considered. 

The  objective  of  this  work  is  to  apply  streamline  simulator 
techniques  to  the  Second  Wall  Creek  Sand  and  predict  the  results  of  a 
waterflood  under  certain  operating  conditions. 


CHAPTER  II 
THE  STREAMLINE  MODEL 

The  theory  and  development  of  the  streamline  model  have  been 
rigorously  described  by  LeBlanc  and  Caudle(2),   Rust(3),    Wessels(4) 
and  Kazmann(5).     A  review  of  the  assumptions  and  the  mathematics 
is  presented  here. 

2.  1  Basic  Assumptions 

The  basic  assumptions  used  in  formulating  the  streamline  model 
are  as  follows : 

1.  Negligible  gravitational  effects  within  the  reservoir.     There 
is  no  gravity  segregation  of  fluids. 

2.  Steady-state  flow  within  the  reservoir.     Incompressible 
reservoir  fluids.     Operations  are  at  nearly  constant  pressure, 

3.  Homogeneous,    isotropic  reservoir. 

4.  Uniform  dispersion  of  free  gas. 

5.  Oilbank  formation. 

The  simulator  used  is  designed  specifically  for  oil  bank  build-up 
and  will  not  accurately  predict  the  results  of  other  operations  such  as 
pressure  maintenance. 


2.  2  Fundamental  Equations 

The  starting  point  for  the  development  of  the  streamline  model  is 
the  continuity  equation.     This  partial  differential  equation  is  derived  from 
a  three  dimensional  material  balance  and  is  written  as: 

d  (pu    )           3  (pu  )          3  (pu  )  an 
a +      — *^ +     *Z -     9    37-  (2.2.1) 
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The  volumetric  flux  may  be  expressed  by  Darcy's  Law  in  the 
differential  form  as: 


_     q  k     d$  to   o   o\ 

u    =    A    =     "    ~    dT  (2'2'2) 


When  the  volumetric  flux  components  in  the  x,   y  and  z  directions,    as 
described  by  Darcy's  Law,   Equation  2.2.2,   are  substituted  into 
Equation  2.2.  1,   the  following  equation  results: 

a     (p\i  a*\        a_    (p_\  aj>\        _L  (p_^±  5*\  =  6  1L 

ax      V  |JL       Bx/  3y      ^       dy/  oz     \  H      oz/  ?     dt  (2.2.3) 

By  using  the  assumptions  listed  previously,   Equation  2.2.3 
may  be  simplified.     The  resulting  equation  may  be  written  in  the  form 
of  a  La  Place  equation. 

JLE     +    !_R     =      0  (2.2.4) 
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The  solution  of  Equation  2.2.4  yields  the  pressure  at  any  point 
in  a  porous  medium  for  which  the  assumptions  made  are  valid.     Knowing 
the  pressure  at  all  points  in  a  porous  medium,   the  average  velocity  of  a 
fluid  particle,    and  thus  its  movement,    may  be  computed. 

By  transforming  to  cylindrical  coordinates  and  allowing  radial 
flow  about  a  point,    Equation  2.  2.4  may  be  solved  for  the  flow  potential 
in  a  horizontal,    isotropic,   homogeneous  porous  medium.     Restoring 
the  resulting  equation  to  rectangular  coordinates  and  using  the  rule  of 
superposition  for  a  multi-well  system,   the  solution  of  Equation  2.2.4 
becomes 


p(x,y)      =      pm       -      ^      ]T       q.    ton         (x-x.)      +    (y-y.) 


i  =  1 


(2.2.5) 


where  p      is  the  mean  reservoir  pressure, 
m 

In  a  porous  medium,   the  relationship  of  flux  to  velocity  is 

-       u       -  k        dP  o    9     A\ 

v  -  j  ~  -  *n  df  (2-2-6) 

By  differentiating  Equation  2.  2.  5  with  respect  to  x  and  substituting  the 

resulting  derivative  into  Equation  2.2.  6  the  velocity  in  the  x  direction 

in  a  rectangular  coordinate  system  becomes: 

n                                x-x. 
Vx(x'y)    =     "    2^h      51      qi      7 N2— 7 y  (2'2-7) 
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In  a  similar  manner,   the  velocity  in  the  y  direction  becomes: 

i         n  y  -y. 

Vy(X'y)     =     "    2^h        5        ^       (  Y     +    (  Y  <2'2'8> 

Using  vector  addition,   the  velocity  of  a  particle  of  fluid  is: 


2  2 
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2.3  Finite  Difference  Equations  for  Movement 

In  order  to  apply  Equations  2.  2.  7  and  2.  2.  8  to  the  streamline 
model,   it  is  assumed  that  the  velocity  remains  constant  during  any  one 
small  distance  increment,   As.     Then  the  time  increment,   At,    is 
calculated  by: 

At      =    —  (2.3.1) 

v  , 

From  the  particles  starting  point  (x    ,   y   ),   it  will  move  to  a  new  point 

s        s 

calculated  as  follows : 

x  =    x      +    v     (x    ,y   )At 

s  + 1  s  x       s      s 

(2.3.2) 

y   j.  i    "   y     +  v    (x  ,  y  )At 

s  + 1  s  y       s      s 

By  starting  on  the  wellbore  of  an  injection  well  and  successively 
applying  Equations  2.2.7  through  2.3.2,    a  trace  of  any  particle's  move- 
ment may  be  made.     This  trace,    called   a  streamline,   is  followed  until  it 
reaches  a  production  well  whereupon  it  is  terminated.     The  travel  time 
for  a  single  streamline  is  the  summation  of  the  time  increments 


corresponding  to  the  number  of  steps  taken  along  the  streamline. 

2.  4         Streamtube  Concept 

In  actual  computation,   a  representative  number  of  streamlines  are 
chosen  to  represent  all  fluids  emanating  from  the  source.     An  example  of 
the  streamlines  for  a    bounded  pattern  are  shown  in  Figure  1.     A  few 
typical  streamlines  are  shown  in  Figure  2.     By  dividing  the  distance 
between  two  adjacent  streamlines  into  two  equal  parts,    the  streamline 
may  be  said  to  represent  the  flow  of  all  the  particles  within  the  dividing 
lines.     A  streamtube  has  thus  been  described  and  is  shown  in  Figure  3(3). 
As  a  fluid  particle  is  traced  along  its  streamline,   the  volume  of  fluid 
produced  at  the  streamline's  production  well  can  be  calculated  and 
accumulated. 

2.5  Multiple  Fluid  Flow 


The  streamline  model  described  so  far  deals  with  a  single  fluid. 
In  a  waterflood  operation,   water  is  injected  into  the  reservoir  to  displace 
the  oil  and  gas.     Therefore  as  many  as  three  fluid  regions  may  be  present 
in  the  reservoir  at  any  time  and  may  each  have  a  different  mobility. 

The  first  assumption  for  the  multiple  fluid  system  is  that  the  oil 
is  miscibly  displaced  by  the  water  in  a  piston-like  manner.     The  oil  is 
banked  up  ahead  of  the  water  such  that  all  oil  displacement  takes  place  at 
the  leading  edge  of  the  water  front.     A  diagram  of  the  concentrations  is 
shown  on  page  9. 
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Figure  1    Streamlines  for  a  Bounded  Pattern. 
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Figure  2     Typical  Streamlines. 
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Figure  3     Typical  Stream  tube 
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With  the  assumption  of  piston-like  displacement  the  displacing  fluid  inter- 
face can  be  easily  followed  along  the  center  pathline  of  a  streamtube. 

The  most  significant  approximation  used  in  extending  the  stream- 
line model  for  multiple  fluids  is  that  the  pathlines  will  be  the  same  as  the 
streamlines  generated  for  a  single  fluid.     This  approximation,   while  not 
correct,   has  been  found  to  be  reasonable  except  at  extreme  values  of 
mobility  ratio(3). 

As  the  fluids  proceed  along  the  streamtube  the  flow  rate  of  the 
fluid  stream  will  change  because  of  the  different  fluid  region  mobilities. 
To  correct  for  the  changing  fluid  flow  rate  the  concept  of  conductivity 
ratio  is  applied.     Conductivity  and  conductivity  ratio  are  developed  and 
explained  by  Caudle (6).     Briefly  the  conductivity  ratio  is  a  comparison 
of  the  conductivity  of  the  medium  at  some  time  when  more  than  one  fluid 
is  flowing  to  the  conductivity  of  the  medium  initially  when  only  one  fluid 
was  present.     The  general  equation  for  the  conductivity  ratio  is 


y  = 


<Pl  -p0) 
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(2.5.1) 
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The  velocity  Equations  2.2.7  and  2.2.8  are  multiplied  by  the  conduc- 
tivity ratio  to  correct  for  multiphase  flow. 

2.  6  Modification  for  Thickness  Variation 

An  empirical  modification  to  allow  for  varying  thickness  has  been 
developed  by  Wessels(4).     This  modification  is  useful  if  the  thickness  at 
each  point  in  the  reservoir  can  be  defined.     A  system  of  grid  squares 
with  assigned  thicknesses  is  used  to  describe  the  reservoir  thickness 
variation. 

The  modifications  for  thickness  variation  are  as  follows: 


Vx  =    "  2TT<j) 


n 


X  -X. 


(x-x.)2   +     (y-y.y 


(2.6.2) 


y  2TT<f> 


n 


Y-Yi 


(x-xj  +     (y-y.y 


(2.6.3) 


The  factor  (h   /h.)  is  an  approximation  factor  which  reflects  an  increasing 

r*        J 

velocity  as  the  thickness  decreases  and  a  decreasing  velocity  as  the 
thickness  increases.     The  diagram  below  shows  the  relationship  between 
reservoir  thickness  and  velocity  for  a  single  streamline. 
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2.7  Mathematical  Bounding  of  the  Streamline  Model 

An  imaging  technique  developed  by  Lin  (7)  which  can  determine 
the  necessary  rates  for  a  preselected  set  of  image  wells  that  will 
confine  a  system  of  streamlines  within  a  given  area  was  used  for  this 
study.     The  technique  uses  a  least  squares  method  to  determine  the 
image  well  flow  rates  so  that  a  no-flow  boundary  is  established  between 
each  of  the  reservoir  bounding  points. 


CHAPTER  III 

APPLICATION  OF  THE  STREAMLINE  MODEL  TO  THE 
SECOND  WALL  CREEK  SAND 


3.  1  Reservoir  History 

President  Wilson's  Executive  Order  of  April  30,    1915  designated 
the  Teapot  Dome  area  in  Wyoming  as  Naval  Petroleum  Reserve  No.    3 
(NPR  3).     On  April  7,    1922,    the  reserve  was  leased  to  the  Mammoth  Oil 
Company  for  the  purpose  of  exploitation.     Following  a  decision  of  the  U.S. 
Supreme  Court  all  of  the  producing  wells  on  the  reserve  were  shut-in  on 
December  31,    1927.     The  field  remained  shut-in  until  1951-1953  when  an 
exploratory  program  was  initiated;    thereafter  it  was  shut-in  until  1958 
when  an  offset  drilling  program  was  instituted  to  protect  against  drainage 
by  adjacent  operators. 

The  wells  are  distributed  unevenly  in  the  Second  Wall  Creek  Sand 
because  development  was  started  at  the  north  end  of  the  reserve  and  was 
stopped  prior  to  completion.     Based  on  approximate  surface  area  the 
northern  third  of  the  reserve  has  69  wells,   the  middle  third   has  23  wells 
and  the  southern  third    has   13  wells. 

Approximately  3,  600,  000  barrels  of  oil  were  produced  prior  to 
December  31,    1927.      1,200,000  barrels  of  oil  have  been  produced  since 
1951  with  most  of  this  production  occuring  after  1958. 
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3.2  Location  of  Oil  and  Gas 

Figure  4  shows  the  approximate  locations  of  the  gas  cap,   oil  zone 
and  water  in  the  Second  Wall  Creek  Sand(8). 

The  Second  Wall  Creek  Sand  just  beyond  the  northwestern  boundary 
of  NPR  3  has  been  under  waterflood  for  several  years.  The  Navy  has  been 
producing  wells  along  this  boundary  to  protect  against  drainage.  Conse- 
quently a  portion  of  the  Naval  Reserve  has  been  swept  by  the  waterflood. 
To  allow  for  this  production,  170  acres  along  the  northwestern  boundary 
were  omitted  from  the  area  which  was  modeled.  The  total  area  modeled 
was  3000  acres. 

3.3  Reservoir  Properties 

Reservoir  data  and  operating  conditions  were  obtained  from 
Tesoro  Petroleum  Corporation  and  U.   S.    Navy  drawings  of  NPR  3 
(8,    10,    11),    Geological  Survey  and  U.   S.    Navy  publications   (1,    9), 
U.   S.   Navy  correspondence  (12,  13)  and  various  core  data  and  logs  pro- 
vided by  the  Officer  in  Charge  of  Naval  Petroleum  Reserve  No.    3. 

Table  I  lists  the  reservoir  and  fluid  properties  used  for  the 
modeling  study.     A  discussion  of  significant  individual  items  follows. 

Porosity  and  Saturations     The  average  porosity  and  connate  water 
saturation  used  are  the  same  as  those  used  in  previous  estimates  of 
waterflooding  at  the  north  end  of  the  reservoir  (12,  13).     No  significant 
difference  was  noted  between  these  figures  and  the  core  data  examined. 
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Figure  4     Location  of  Oil  and  Gas  in  the  Second  Wall  Creek  Sand. 
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TABLE  I 
RESERVOIR  AND  FLUID  PROPERTIES 

Porosity  (Avg),   %  17.5 

Connate  Water  Saturation,   %  35 

Gas  Saturation,   Initial,   %  5 

Gas  Saturation,   Residual,   %  0 

Oil  Saturation,    Initial,   %  60 

Oil  Saturation,   Residual,   %  20 

Major  Fault  Blocks  10 

Net  Thickness,   Ft.  21-60 

Reservoir  Area,   Acres  3000 

Depth  (Avg),   Ft.  2800 

Temperature,   °F  140 

Pressure,   PSI  1100 

Permeability  (Avg),  MD  10 

Formation  Volume  Factor,   Oil  1.2 

Formation  Volume  Factor,   Water  1.0 

Oil  Gravity,    °API  3  8 

Oil  Viscosity,   CP  6 

Water  Viscosity,  CP  .7 

Gas  Viscosity,   CP  .  12 

Water/ Oil  Mobility  Ratio  4 

Gas  /Oil  Mobility  Ratio  65 

Gas  Oil  Ratio,     SCF/BBL  1000 

Solution  Gas  Oil  Ratio,   SCF/BBL  50 
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For  this  study  it  was  assumed  that  the  waterflood  would  start  at 
the  very  beginning  of  the  production.     The  model  was  run  with  an  initial 
gas  saturation  of  5%  to  allow  oil  bank  build-up  to  occur  and  to  show  the 
effect  of  fill-up  time.     The  initial  oil  saturation  used  was  60%. 

The  residual  gas  saturation  was  set  at  zero  because  it  was 
assumed  that  all  of  the  free  gas  would  either  be  displaced  or  put  back 
into  solution  as  the  oil  bank  passed. 

The  residual  oil  saturation  after  the  waterflood  passed  was  set 
at  20%.     This  is  almost  twice  the   12%  average  residual  oil  saturation 
found  in  the  cores  available.     The  residual  oil  saturation  in  the  cores 
was  considered  an  absolute  minimum  which  would  be  difficult  to  achieve 
in  the  reservoir. 

Fault  Blocks       The  Second  Wall  Creek  Sand  is  heavily  faulted  and 
was,   for  the  presentation  of  well  productivity  data,    divided  into  about 
twenty  possible  fault  blocks(12).     In  order  to  model  the  field,    taking  into 
account  the  faulting,   the  field  was  divided  into  ten  sections  along  the 
largest  faults.     It  was  assumed  that  the  faults  were  sealed  and  the 
sections  were  then  modeled  independently.     This  had  the  added  advantage 
of  making  the  modeled  reservoirs  smaller  and  easier  to  manage.     Also, 
properties  which  vary  within  the  reservoir  could  be  easily  changed  from 
section  to  section.     Figure  5  shows  the  configuration  of  the  ten  fault 
blocks  chosen. 
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Boundary  Naval  Petroleum 
Reserve  No.    3 


--    Oil/Water  Contact 


Scale  :1"=  4000' 


Figure  5    Fault  Blocks  Used  in  Modeling  the  Second  Wall  Creek  Sand. 
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Net  Sand  Thickness     Thin  beds  of  shale  are  common  throughout 
the  Second  Wall  Creek  Sand,   but  they  thicken  toward  the  south  and  east 
and  split  the  sand  into  several  distinct  layers.     Some  wells  have  three  or 
more  layers  of  sand  separated  by  shale  beds  from  two  to  twenty -five 
feet  thick(12)„ 

Combining  the  information  available  from  the  cores  with  the 
general  trends  mentioned  above,    a  percentage  discount  factor  was  used 
for  each  fault  block  to  reduce  the  gross  sand  thickness  to  a  net  sand 
thickness.     Table  II  shows  the  discount  factors  that  were  applied. 
Figure  6  shows  the  gross  thickness  of  the  Second  Wall  Creek  Sand(lO). 

Permeability    The  temperature,   pressure  and  permeability 
were  used  to  determine  the  fluid  and  fluid  flow   properties.     The  perme- 
ability itself  is  not  used  in  the  streamline  model  calculations  because  the 
calculations  are  made  at  a  constant  flow  rate.     This  means  that  the 
injection  pressure  will  have  to  be  great  enough  to  make  the  injectors 
flow  at  the  required  rate.     Rate  selection  will  be  discussed  in  the 
next  chapter. 

Gas  -Oil  Ratios       Gas -oil  ratios  vary  greatly  depending  on  where 
in  the  field  the  well  is  located.     The  gas -oil  ratio  and  solution  gas -oil 
ratio  are  used  by  the  computer  program  to  indicate  the  arrival  of  the 
oil  bank.     Representative  values  of  the  two  ratios  were  used  but  they  do 
not  represent  a  rigid  material  balance. 
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TABLE  II 

DISCOUNT  FACTORS  APPLIED  TO 
GROSS  SAND  THICKNESS 
TO   OBTAIN  NET  SAND 
THICKNESS 


Fault  Block  Factor 

1  0,8 

2  0,7 

3  0,7 

4  0,8 

5  1.0 

6  1.0 

7  0.8 

8  0.7 

9  0,  66 
10  0o66 
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Section  34 
Fault 


■"•  Boundary  Naval  Petroleum 
Reserve  No.    3 

--Oil/Water  Contact 
Scale  :1"=  4000' 


Figure  6   Isopachous  Map  of  the  Second  Wall  Creek  Sand 
Gross  Sand  Thickness. 
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3.  4  Operating  Conditions 

Table  III  is  a  listing  of  the  operating  conditions  used  for  the 
simulation. 

TABLE  III 

OPERATING    CONDITIONS  FOR  THE  WATERFLOODING 
OF  THE  SECOND  WALL  CREEK  SAND 


No.    of  Production  Wells 

135 

No.    of  Injection  Wells 

132 

Well  Pattern 

5  -Spot/  Line  Drive 

Well  Spacing,   Acres 

10 

Maximum  Oil  Production  Rate, 

STB/Well/Day 

46-300 

Oil  Production  Rates  During  Fill  Up, 

STB/Well/Day 

3-25 

Total  Water  Injection  Rate,   STB/Day 

23,000 

Field  Development    Based  on  a  proposed  Navy  Development 
Plan(ll)  the  field  was  assumed  to  be  drilled  on  ten  acre  spacing. 

Well  water  injectivity  was  assumed  to  be  in  the  same  range  as 
the  well's  initial  oil  productivity.     Therefore,    approximately  a  one-to- 
one  ratio  of  injectors  to  producers  was  maintained. 

The  injection-production  well  pattern  used  was  a  combination 
five -spot  and  line  drive.     The  injection  wells  were  generally  placed 
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along  the  outside,    nearest  the  oil/water  contact,   in  a  line  drive  pattern  to 
force  the  oil  toward  the  center.     Then,   where  area  allowed,    injection  wells 
were  placed  within   the  producers  in  a  five -spot  pattern. 

Production  Well  Rates     In  the  streamline  model  the  oil  production 
rate  of  a  well  during  the  fill-up  time  is  set  equal  to  the  present  rate  to 
which  the  well's  oil  production  has  declined.     After  oil  bank  breakthrough 
the  oil  production  rate  of  a  well  is  set  equal  to  the  initial  or  maximum 
oil  rate  of  that  well. 

The  initial  rates  for  each  production  well  were  determined  from 
a  composite  drawing  of  the  Second  Wall  Creek  Sand  production(8).     Since 
the  reservoir  will  be  waterflooded  from  the  start,    the  oil  rates  during 
fill -up  were  set  at  10%  of  each  well's  initial  rate.     At  reservoir  conditions, 
the  total  water  injection  rate  was  set  equal  to  the  total  oil  production  rate 
after  fill-up  to  satisfy  the  steady  state  assumption. 

The  individual  injection  rates  were  assigned  to  the  injection  wells 
in  such  a  manner  to  contain  the  streamlines  as  much  as  possible  to  the 
oil  bearing  area.     Some  flow  into  the  gas  cap  and  the  watertable  cannot 
be  stopped  but  proper  balancing  minimized  the  losses. 

Figures  7  and  8  show  the  well  pattern  used  and  the  streamlines 
generated  by  the  model  reservoir.  See  Appendix  A  for  a  listing  of  the 
individual  well  data  used. 
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X   Injection  Wells 
O   Production  Wells 
Scale  :1"^  2300' 


Section  34 
Fault 


Figure  7     Model  Generated  Streamlines -Northern  Half   of 
Second  Wall  Creek  Sand. 
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Section  34 
Fault 


X    Injection  Wells 
O    Production  Wells 
Scale  :1"=  2300' 


Figure  8     Model  Generated  Streamlines -Southern  Half  of 
Second  Wall  Creek  Sand. 


CHAPTER  IV 
RESULTS  OF  THE  STREAMLINE  MODEL  SIMULATION 

4.  1  Predicted  Production  History- 

Figures  9,    10  and  11  show  the  predicted  response  of  the  Second 
Wall  Creek  Sand  to  a  waterflood  operated  as  previously  described.     It  is 
predicted  that  18,  000,  000  barrels  of  oil  will  be  recovered  in  five  years 
and  30,  000,000  barrels  of  oil  after  ten  years.     "Water  breakthrough  occurs 
around  1,  000  days  and  the  water -oil  ratio  climbs  to  over  five  after  ten 
years.     Ultimate  recovery  at  a  water-oil  ratio  of  twenty  is  predicted 
to  be  around  32,  000, 000  barrels  and  takes  twelve  years. 

4.2  Effects  of  Major  Assumptions  on  Recovery  Predictions 

Although  the  oil  recovery  predicted  agrees  with  the  present 

reserve  estimate  of  28,  800,  000  barrels (9),   the  predictions  of  both 

quantity  and  time  are  only  as  good  as  the  assumptions  made. 

In  order  to  provide  a  proper  perspective,   the  variables  and 

assumptions  which  could  most  effect  the  results  will  be  discussed. 

Heterogeneity     The  model  was  run  assuming  that  the  reservoir 
sand  in  each  fault  block  was  homogeneous.     Frequently,   however,   hori- 
zontally bedded  sandstones  exhibit  different  permeabilities  among  the 
individual  stratum.     The  presence  of  high  permeability  strata  will 
reduce  the  efficiency  of  the  waterflood.     The  high  permeability  zones 
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water-out  faster  than  the  rest  of  the  reservoir  and  then  increase  the 
water  oil  ratio  much  sooner  than  expected.     If  high  permeability  strata 
are  shown  to  be  a  significant  factor  the  predictions  presented  here  will 
be  optimistic. 

Initial  Gas  and  Oil  Saturation     The  presence  of  free  gas  in  the 
oil  zone  of  the  reservoir  means  that  there  will  be  a  certain  time  lag, 
fill -up  time,   before  the  oil  bank  will  reach  the  production  wells.     This 
time  lag  is  approximately  equal  to  the  time  required  to  inject  enough 
water  to  fill  the  gas  space  and  the  space  of  the  oil  that  is  produced.     The 
initial  gas  saturation  chosen  to  represent  the  reservoir  was  5%.     If  the 
quantity  of  gas  present  in  the  reservoir  is  actually  less  than  5%  the  oil 
production  rate  will  increase  earlier  and  if  the  space  is  assumed  to 
be  filled  with  oil,   more  oil  will  be  recovered.     If  the  gas  saturation  is 
greater  than  5%  the  opposite  will  be  true. 

Production  and  Injection  Rates      The  effect  of  increasing  or 
decreasing  the  injection  and  production  rates  would  be  to  condense  or 
expand  the  time  scale.     If  for  example  it  was  found  that  the  reservoir 
could  operate  at  injection  and  production  rates  50%  greater  than  those 
assumed  here,   the  recovery  that  would  have  taken  1,  500  days  could  be 
achieved  in  1,000  days.     See  Appendix  B  for  an  example  of  rate  adjust- 
ment to  obtain  uniform  water  breakthrough  times. 


CHAPTER  V 
CONCLUSIONS  AND  RECOMMENDATIONS 


5.  1  Conclusions 


The  waterflood  of  the  Second  Wall  Creek  Sand  will  increase  the 
total  oil  production  to  around  32,  000,  000  barrels.     This  is   15,  000,  000 
barrels  over  that  which  could  be  expected  from  solution  gas  drive  primary- 
production.     At  the  operating  rates  used  the  waterflood  will  require  around 
twelve  years. 

5.2  Recommendations 

1.  It  is  recommended  that  the  waterflooding  of  the  Second  Wall 
Creek  Sand  be  planned  and  cost  evaluations  made.     The  streamline  model 
developed  here  is  well  suited  for  studying  development  and  operating 
considerations  so  that  project  economics  can  be  evaluated. 

2.  Since  the  effects  of  high  permeability  strata  on  secondary 
recovery  operations  are  always  adverse  and  sometimes  severe,    it  is 
recommended  that  the  possible  existence  of  such  strata  be  thoroughly 
investigated  during  waterflood  planning  and  development. 
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NOMENCLATURE 


Symbol 

Definition 

A 

Cross  Sectional  Area 

h 

Reservoir  thickness 

k 

Permeability 

1 

Length;  number  of  fluids 

n 

Number  of  wells 

P 

Pressure 

q 

Flow  rate 

As 

Distance  increment  alon 

s  Saturation 

At  Time  step  along  path  line 

t  Time 

u  Flux 

v  Average  fluid  particle  velocity 

x,y,  z  Rectangular  coordinates 


y 

Conductivity  ratio 

X 

Fluid  mobility 

n 

Viscosity 

p 

Density 

o 

Porosity 

$ 

Potential 

31 


32 


Subsc 

ripts 

cw 

Connate  water 

gi 

Initial  gas 

gr 

Residual  gas 

i 

Index  of  wells;  in 

j  Index  of  particle  location 

k  Index  of  fluids 

m  Mean 

o  Oil;  out 

oi  Initial  oil 

or  Residual  oil 

s  Starting  point  of  fluid  particle 

p  Production 

w  Water 

x  In  the  x-direction 

y  In  the  y-direction 

z  In  the  z  -direction 


APPENDICES 


APPENDIX  A 
INDIVIDUAL  WELL  DATA 
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TABLE  IV 
WELL  DATA— FAULT  BLOCK  1 
Refer  to  Figure  12 


Well 

Rate, 

, STB/Day 

Thickness 

Well 

Rate 

,  STB 

/Day 

Thickness 

No. 

(+In; 

i,  -Prod) 
87. 

Ft. 

No. 

(+Irv 

j,  -Prod) 
96 

Ft. 

1 

29 

21 

32 

2 

54. 

27 

22 

-96 

32 

3 

-90. 

30 

23 

48 

24 

4 

52. 

26 

24 

-72 

24 

5 

-87. 

29 

25 

-72 

24 

6 

93 

31 

26 

-78 

26 

7 

56 

28 

27 

84 

28 

8 

-90 

30 

28 

-90 

30 

9 

-96 

32 

29 

96 

32 

10 

78 

26 

30 

-96 

32 

11 

-81 

27 

31 

-96 

32 

12 

90 

30 

32 

48 

24 

13 

-96 

32 

33 

-72 

24 

14 

96 

32 

34 

72 

24 

15 

48 

24 

35 

-72 

24 

16 

72 

24 

36 

78 

26 

17 

-78 

26 

37 

-84 

28 

18 

-81 

27 

38 

90 

30 

19 

87 

29 

39 

-96 

32 

20 

-93 

31 

40 

96 

32 
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TABLE  IV  (Continued) 


Well 

Rate, 

STB/Day 

Thickness 

Well 

Rate 

,  STB/Day 

Thickness 

No. 

(+Inj 

,  -Prod) 
48 

Ft. 

No. 
61 

(+Inj,  -Pi 
-90 

rod) 

Ft. 

41 

24 

30 

42 

-72 

24 

62 

81 

27 

43 

72 

24 

63 

48 

24 

44 

-72 

24 

64 

36 

24 

45 

72 

24 

65 

72 

24 

46 

-78 

26 

66 

-78 

26 

47 

84 

28 

67 

93 

31 

48 

-90 

30 

68 

-90 

30 

49 

96 

32 

69 

84. 

28 

50 

-96 

32 

70 

-84 

28 

51 

96 

32 

71 

48. 

24 

52 

48 

24 

72 

48. 

24 

53 

-72 

24 

73 

-78 

26 

54 

-72 

24 

74 

-78 

26 

55 

-72 

24 

75 

36 

24 

56 

90 

30 

76 

40 

24 

57 

-90 

30 

77 

48 

24 

58 

90 

30 

59 

-90 

30 

60 

90 

30 

Production  well  oil  rates  during  fill-up  equal  10  STB/Day 
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Figure  12  Well  Locations  for   Fault  Blocks  1    and  2. 
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TABLE   V 
WELL  DATA— FAULT  BLOCK  2 
Refer  to  Figure  12 


Well 

Rate,  STB 

/Day 

Thickness 

Well 

Rate,  STB /Day 

Thickness 

No. 

(+Inj,  -Prod) 

Ft. 

No. 
20 

(+Inj,  -Prod) 

'1- 

Ft. 

1 

28 

29.4 

2 

* 

30.  1 

21 

30.  8 

3 

* 

31.5 

22 

-175 

29.4 

4 

32.2 

23 

-180 

30.  1 

5 

-175 

28.7 

24 

28.7 

6 

-180 

30.  1 

25 

29.4 

7 

-195 

32.2 

26 

-155 

25.9 

8 

-205. 

34.3 

27 

-175 

28.7 

9 

-175 

28.7 

28 

* 

24.5 

10 

-185 

30.8 

29 

-135 

22.4 

11 

-200 

32.9 

30 

-147 

24.  5 

12 

-210 

35. 

31 

* 

25.  9 

13 

•A. 

29.4 

32 

-160 

26.6 

14 

*>* 
"1* 

30.8 

33 

V 

28.7 

15 

* 

33.6 

34 

-175 

29.4 

16 

* 

35. 

35 

* 

21.0 

17 

-175 

29.4 

36 

23.  1 

18 

-185 

30.8 

37 

:'; 

25.  9 

19 

-200 

32.9 

38 

28.7 

Production  well  oil  rates  during  fill-up  equal  20  STB/Day. 

'"Program  calculates  injection  well  rates  based  on  reservoir  thickness 
and  steady  state  material  balance. 
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TABLE  VI 
WELL  DATA—FAULT  BLOCK  3 
Refer  to  Figure  13 


Well 

Rate,  STB /Day 

Thickness 

No. 

(+Inj,  -Prod) 
140 

Ft. 

1 

35. 

2 

140 

35 

3 

140 

35 

4 

-178.5 

35.7 

5 

-182 

36.4 

6 

-178.5 

35.7 

7 

218.4 

36.4 

8 

226,2 

37.7 

9 

226.2 

37.7 

10 

-182 

36.4 

11 

-192. 5 

38.5 

12 

-203, 

40.6 

13 

-188.5 

37.7 

14 

175 

35. 

15 

182 

36.4 

16 

203 

40.6 

17 

-210 

42. 

18 

-182 

36.4 

Production  well  oil  rates  during  fill-up  equal  20  STB/Day. 
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Figure  13     Well  Locations  for  Fault  Blocks   3-5. 
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TABLE  VII 

WELL  DATA—FAULT  BLOCK  4 

Refer  to  Figure   13 


Well 

Rate, 

STB/Day 

Thickness 

No. 

(+Inj 

,  -Prod) 
88 

Ft. 

1 

44 

2 

-100 

50 

3 

90 

45 

4 

-100 

50 

5 

86 

43 

6 

-94 

47 

7 

48 

48 

8 

-48 

48 

9 

-48 

48 

10 

48 

48 

Production  well  oil  rates  during  fill-up  equal  10%  of  initial  rate 
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TABLE  VIII 

WELL  DATA— FAULT  BLOCK  5 

Refer  to  Figure  13 

Well         Rate,  STB/Day         Thickness 
No.  (-Unj,  -Prod)  Ft. 

1  -50  50 

2  *  50 

3  -50  50 

4  *  48 

5  -50  50 

6  *  50 

7  -46  46 

8  *  46 

Production  well  oil  rates  during  fill-up  equal  3  STB/Day,, 

*  Program  calculates  injection  well  rates  based  on  reservoir  thickness 
and  steady  state  material  balance. 
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TABLE  IX 
WELL  DATA— FAULT  BLOCK  6 
Refer  to  Figure  14 


Well 

Rate,  STB/ 

'Day 

Thickness 

No. 

(+Inj,  -Prod) 
260 

Ft. 

1 

56 

2 

240 

52 

3 

245 

53 

4 

245 

53 

5 

225 

49 

6 

-340 

58 

7 

-300 

60 

8 

-300 

60 

9 

-300 

60 

10 

-285 

57 

11 

-255 

51 

12 

-275 

55 

13 

300 

60 

14 

-300 

60 

15 

-270 

54 

16 

250 

48 

17 

305 

52 

18 

-275 

55 

Production  well  oil  rates  during  fill-up  equal  25  STB/Day. 
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FAULT  BLOCK    6 


Section    3 


FAULT  BLOCK   7 


Section   4 


9 


FAULT  BLOCK  8 


Section    10 


x    Injection  Well 
O    Production  Well 


Figure  14     Well  Locations  for   Fault  Blocks  6-8 


• 


. 
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TABLE  X 
WELL  DATA— FAULT  BLOCK  7 
Refer  to  Figure  14 


Well 

Rate 

,STB 

/Day 

Thickness 

Well 

Rate,  STB/Day 

Thickness 

No. 

(+Inj,  -P: 

234 

rod) 

Ft. 

No. 
13 

(+Inj,  -Pi 
-282 

rod) 

Ft. 

1 

39 

47 

2 

-264 

44 

14 

258 

43 

3 

234 

39 

15 

-288 

48 

4 

-258 

43 

16 

276 

46 

5 

234 

39 

17 

-288 

48 

6 

-258 

43 

18 

270 

45 

7 

228 

38 

19 

-282 

47 

8 

-240 

40 

20 

264 

44 

9 

222 

37 

21 

-276 

46 

10 

-228 

38 

22 

258 

43 

11 

-270 

45 

23 

-270 

45 

12 

252 

42 

24 

258 

43 

Production  well  rates  during  fill-up  equal  25  STB/Day. 
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TABLE  XI 

WELL  DATA— FAULT  BLOCK  8 

Refer  to  Figure  14 


Well 

Rate, 

STB/Day 

Thickness 

No. 

(+Inj 

i,  -Prod) 
124 

Ft. 

1 

31 

2 

-124 

31 

3 

124 

30 

4 

-124 

31 

5 

-168 

42 

6 

168 

42 

7 

-168 

42 

8 

160 

40 

9 

-168 

42 

10 

156 

39 

Production  well  oil  rates  during  fill-up     equal  15  STB/Day. 
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TABLE  XII 
WELL  DATA— FAULT  BLOCK  9 
Refer  to  Figure   15 


Well 

Rate , , 

STB/Day 

Thickness 

Well 

Rate,, 

STB/Day 

Thickness 

No. 

(+Inj 

,  -Prod) 

Ft. 

No. 
14 

(+Inj 

,  -Prod) 
-150 

Ft. 

1 

27.7 

40 

2 

# 

30.4 

15 

* 

40 

3 

-120 

29. 

16 

-150 

40 

4 

-160 

31.7 

17 

-200 

40 

5 

1* 

30.4 

18 

* 

40 

6 

-I- 

35.6 

19 

-150 

40 

7 

-135 

33. 

20 

V 

40 

8 

-200 

40. 

21 

-150 

40 

9 

•A. 

35.6 

22 

-150 

40 

10 

40 

23 

40 

11 

-155 

37.6 

24 

-1,. 

40 

12 

-200 

40. 

25 

-150 

40 

13 

i'z 

40 

Production  well  oil  rates  during  fill -up  equal  15  STB/Day. 

*  Program  calculates  injection  well  rates  based  on  reservoir  thickness 
and  steady  state  material  balance. 


48 


FAULT  BLOCK    9 


Section    15 


FAULT  BLOCK    10 


i       2^      AP 
15  \„o  *f 


11 


Section    14 


X    Injection  Well 
O    Production  Well 


Figure  15    Well  Locations  for  Fault  Blocks  9  and  10. 
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TABLE  XIII 
WELL  DATA— FAULT  BLOCK  10 
Refer  to  Figure  15 


Well 

Rate 

,  STB/Day 

Thickness 

Well 

Rate 

,  STB 

/Day 

Thickness 

No. 

(+Inj,-Prod) 
-76 

Ft. 

No. 
21 

(+In; 

I,  -Prod) 

Ft. 

1 

38 

40 

2 

* 

35 

22 

-80 

40 

3 

* 

40 

23 

-80 

40 

4 

-76 

38 

24 

-80 

40 

5 

-70 

35 

25 

40 

6 

-80 

40 

26 

-76 

38 

7 

-74 

37 

27 

35 

8 

35 

28 

J, 

'1- 

40 

9 

a* 

40 

29 

-80 

40 

10 

-80 

40 

30 

-78 

39 

11 

-74 

37 

31 

-76 

38 

12 

* 

34 

32 

-72 

36 

13 

1- 

40 

33 

-70 

35 

14 

-150 

40 

34 

'I- 

38 

15 

-150 

40 

35 

-76 

38 

16 

* 

40 

36 

•>* 

37 

17 

-80 

40 

37 

* 

36 

18 

-76 

38 

38 

* 

33 

19 

-72 

36 

39 

-66 

33 

20 

* 

34 

Production  well  oil  rates  during  fill-up  equal  8  STB/Day. 

*  Program  calculates  injection  well  rates  based  on  reservoir  thickness 
and  steady  state  material  balance. 


APPENDIX  B 
RATE  ADJUSTMENT 
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ADJUSTMENTS  OF  INJECTION  AND  PRODUCTION  RATES  TO 
ACHIEVE  UNIFORM  WATER  BREAKTHROUGH 
IN  ALL  FAULT  BLOCKS 

When  the  model  response  was  inspected  section  by  section  it  was 
found  that  some  of  the  sections  were  watering-out  faster  than  others. 
This  situation  is  also  possible  in  the  field.     If  the  injection  and  production 
rates  cannot  be  sufficiently  adjusted,    certain  sections  would  be  watered  - 
out  and  shut-in  while  others  would  be  run  longer  until  ultimate  recovery. 

If  the  rates  can  be  adjusted  enough  the  field  can  be  operated  so 
that  all  of  the  sections  water-out  at  about  the  same  time.     The  well  rates 
used  in  the  model  were  adjusted,   with  the  total  field  injection  rate  held 
constant,   to  determine  the  effect  of  uniform  water  breakthrough  from 
section  to  section. 

The  results,    as  shown  in  Figure  16,   indicate  that  for  the  same 
period  of  operation  the  oil  production  was  increased  and  the  water 
production  was  decreased. 
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APPENDIX  C 

COMPUTER  PROGRAM 

1.  Program  Listing 

2.  Sample  Data  Deck 
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C 
C 
C 
C 
C 

c 
t 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


PROGRAM  OLBNKVTdNPUT, OUTPUT) 


A  STREAMLINE  MOOEL  FOR  WATERFLOOD  SIMULATION 


THIS  MOOEL  IS  SPECIFICALLY  DESIGNED  TO  SIMULATE  A 
"ATERFLOOO  OPERATION  IN  WhICH  An  OIL  BANK  IS  FORMED  AND 
MOVES  THROUGH  THE:  RESERVOIR.   THE  MODEL  15  CAPABLE  OF 
HANDLING  RESERVOIRS  WITH  THICKNESS  VARIATIONS  ANO' 
IRREGULAR  BOUNDARIES. 

REFERENCES: 

LEBLANCtJ.L.*  AND  CAUDLEtB.H.,  A  STREAMLINE  MODEL  FOR 
SECONDARY  RECOVERY*   SOC.  P£T.  EngR.  J.  (I97l)»ll#7. 

RUST.C.B.t  A  STREAMLINE  MOOEL  FOR  OIL  BftNK  BUUOJP  IN 
A  WATERFLOOD*  M.S.  THESIS.  THE  UNIVERSITY  OF  TEXAS  AT 
AUSTIN  (MAY.  1972). 

WESSELSt J.W..  APPLICATION  OF  THE  STREAMLINE  RESERVOIR 
MODEL  TO  RESERVOIRS  HAVING  VARIATIONS  IN  FORMATION 
THICKNESS*  M.S.  THESIS.  THE  UNIVERSI1Y  OF  TEXAS  AT  AUSTIN 
(MAY.  1973). 

DOCUMENTED  BYt  K.E.GOLTZ  (MARCH, 1975) 


*   SCHEMATIC  PROGRAM  FLO*  CHART   • 


I  read  input  data  i 
i 


I   COMPUTE  STREAMLINE  STARTIN3  POINTS.    I 
I    SOURCE  AND  SINK  PRESSURES  I 


I   ITERATE  EACH  TIME  PERIOD    I 


■I   ITERATE  EACH  PRODUCTION  *ELL    I 


I 


'-——I   ITERATE  EACH  STREAMLINE     I 

"l 

-I   COMPUTE  VELOCITIES  AND  TIME  INCREMENTS    I 
I     AT  EACH  POINT  ALONG  STREAMLINE  I 


I 


I   MOVE  FLUID  PARTICLE  ALONG  STREAMLINE  USING   I 
I   FINITE  DIFFERENCE  ApPROX.  OF  DARCYS  LAW      I 
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C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

e 

c 

c 

c 
c 
c 
c 
c 


I 


I   ACCUMULATE  STREAMLINE  PRODUCTION    I 


I 


I  TEST  FOR  BREAKTHROUGH  OF  FLUID 

I  PARTICLES  INTO  PRODUCTION  WELLSt 

I  ADJUST  STREAMLINE  PRODUCTION 
I    RATES  AS  REQUIRED 


I 


NO 


NO 


NO 


-I   TEST  FOR  END  OF  TIME  PERIOD     I 

I  YES 
—  I   TEST  FOR  ALL  STREAMLINES      I 

I  YES 
■I    TEST  FOR  ALL  PRODUCTION  WELLS     I 

I   YES 

I  print/plot  well  and  fie'.o  results  i 
"l 

I   TEST  FOR  TOTAL  TIME   I 
"i   YES 
I   STOP   I 


NO 


»••«••«»•»•»»«•»«•*•»»••««»»•««»••«••««»»»««•»••»•»»•»»»»••»»•« 


INPUT  DATA  DEFINITIONS 


POR  » 
SC*  « 
SOI  i 
SOR  » 
SGR  : 
SGOR 


porosity  (decimal) 
saturation  connate  wat 
saturation  oil  initial 
saturation  oil  residua 
saturation  gas  residua 
:  solution  gas-oil  rati 


NP  s  NUMBER  OF  PRODUCTION  WE 
NI  «  NUMBER  OF  INJECTION  WEL 
NBW  »  NUMBER  OF  IMAGE  WELLS 
NBP  »  NUMBER  OF  POINTS  USED 
NST  s  NUMBER  OF  STREAMLINES 
NXH  e  NUMBER  OF  X-DIRECTlON 
NYH  a  NUMBER  OF  Y-OJRECTION 
NEDIT  «  0  OR  1 

IF  6  OR  BLANK-PROGRAM  WILL 

IF  1   PROGRAM  INTERNAL  BAl 

H  «  AVERAGE  RESRVOIR  THICKNE 

RI  s  RADIUS  OF  STARTING  AND 

DELT  «  TIME  INCREMENT  (DAyS) 

CORRECTION  FACTOR.  SI 


ER  (DECIMAL) 

(DECIMAL) 
L  (DECIMAL) 
L  (DECIMAL) 

0  (SCF/BBL  AT  RES.  CON3HTIONS) 
LLS 
LS 

JSED  TO  BOUNO  PATTERN 
TO  DEFINE  BOUNDARY 
?OR  LOWEST  RATE  PRODUCTION  WELL 
ENTRIES  IN  THICKNESS  GRID 
ENTRIES  IN  THICKNESS  GRID 

BALANCE  PROD/INJ  RATES 
ANCE  FEATURE  IS  DEFEATED 
SS.FT  (USED  FOR  IMAGE  WELLS) 
ENDING  ClRCLESi^T 

FOR  UPDATING  CONSTANT  ORESSJRE 
NCE  DELT  MUST  START  OUT  SMALL  BUT 


62 

63 

6* 

63' 

66 

67 

69 

69 

75 

71 

72. 

73 

7* 

75' 

76 

77 

79 

7? 

8D 

91 

92 

93 

8V 

9  5' 

95 

97 

99 

89 

90 

91 

92 

93 

9ft 

95' 

95 

97 

99 

99 

100 

101 

10> 

103 

10ft 

105' 

105' 

107 

109 

109 

110 

Hi 

112 

113 

ML 

115' 

116 

117 

119 

119 

120 

121 

122 

123 
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C  BECOMES  LESS  CRITICAL  AS  ThE  SOLUTION  PROGRESSES, 

C  THE  VALUE  O'l  DELT  IS  MULTlPlED  BY  TEM  AFTER  THE 

C  OIL  BANK  BREAKS  THROUGH  IN  THE  FIRST  STREAMLINE* 

C  DELP  n  TIME  INCREMENT  FOR  PRINT  AND  PLOT  OAYS) 

C  FIRST  PRINT/PLOT  Is  MADE  AT  DELT.THEN  EVERY  DELP< 

C  TMX  s  MAXIMUMTIME  FOR  SOLUTION  (DAYS) 

C  QMNP  s  LOWEST  PROOJCTION  WELL  RATE  (STB/DAY) 

C  QMNI  a  LOWEST  INJECTION  WELL  RATE  (ST3/DAY/FT) 

C  XMX  s  ABSOLUTE  X  OIMENSION  AT  WHICH  STREAMLINES  AREl  CUT  OFF' 

C  YMX  b  ABSOLUTE  Y  DIMENSION  AT  WHICH  STREAMLINES  AR£l  CUT  OFF. 

C  GM8»OM9»WM8  >  GASfOIL  AND  WATER  MOBILITIES 

C  XW(I),YW(I)  *    X»Y  COORDINATES  OF  WELL  LOCATION 

C  HKI)  ■■  THICKNESS  OF  RESERVOIR  AT  WELL, FT 

C  Q(D  ■  WELLINlTIALi  OIL  PROOJCTION  RATE  (ST3/OAY) 

C  QOU>  «  WELL  PRESENT  (FILL-UP)  OIL  PRODUCTION  RATE  (STB/DAY) 

C  GOR(I)  *  WELL  GAS-OIL  RATIO  (SCF/98L  AT  R£S.  CONDITIONS) 

C  XB(I),YB(I)  «  X»Y  COORDINATES  OF  BOUNDARY  POINTS 

C  XH(I),YH(I)  c  X,Y  COORDINATE  VALUES  FOR  THICKNESS  GQID 

C  HT(I.J)  =  VALUES  OFi  THICKNESS  CORRESPONDING  TO  COORDINATES 

C  OF  THICKNESS  GRID 

C 

C  •••»••♦••#••••»»*•»»•  »»»«»«  •»«»»••••««««  •»•»«»  »»♦»»»»»•»*»*♦»  *>♦»»»» 

c 

COMMON/A A A/  XW ( 125) . Yw < 125) , Hi ( 125) ,Q ( 125) , XB ( 125) , 
IYBU25)  »N8P,NP,NI,NBW,H 

COMMON/BBB/  C( 55, 55) ,0(55) ,01(55) 

COMMON/CCC/  X(8O,20) »Y(80,2o) »X0( 80,20 ) »Yu(30«20>  t 
1WQ(80,20) ,SCL,XMX,YMX,NST,QMNP,RGO<80»20)« 
2SG0R,NPIW,TMEI 

OIMENSION  TSUM(8o»20> • PP < 80. 20) ,Pl (80*20> • CO (80* 20) ♦ 
1CW (80, 20) ,QO(125) ,30R(125) ,00(80*20) » 
2XH(20) ,YH(20) ,HT(20,20) 

READ  530,  P0R,SCW,S0I,S0R*S6RfSG0R 

READ  540,  NP,NI,N8W,NBP,NST 

READ  5*0,  NXH,NYH,NEOIT 

READ  530,  H,RI,OELT,DELP»TMX,QMNP 

READ    530,    QMN|I*XMX,YMX,GMB»OHB,WMB 

DELTT=OELTftlO. 

NCOUNT«0 

PLTFT=0»9»XMX/ YMX* 0.45 

TOTFT»0.7*PLTFT 

NPlWsNP*Nl 

NPIBWzNPlW«NBW 

READ  550,  (XW(I),YW(I),HI(I),Q(I)  tQO(I),GOR(D,Isl,NPIBW) 

READ  560,  (XB(I) ,Y3(l),lsl,N3P) 
C 

C        CONVERT  SURFACE  «ATES  To  RESERVOIR  RATES 
C  FVF'  ■  lt2 

C 

DO  9  I«1,NPIW 

Q(I)»Q(I)*>1.2 
9    QO(I)«30(I)fti.2 

OMNP«QMNP*1.2 

QMMl«OMNI»iig 

PRINT  590 

PRINT  480*  NP,NI,N3W»NBP,NST 

PRINT  490,  H,DELT,TMX,RI 

PRINT  500»  DELPiXMX, YMX, QMNP 

PRINT  510*  P0R,SCW,S0I,S0R,S6R 

PRINT  520*  OMB.WMd.GMB 

PRINT  550t  (XW(I) tYW(I) ,HI(I),Q(I)*QO(I)tG09(I) tI*l,NPIBN) 

PRINT  530*  (XB(I) *YB(I)*Isl,NBP) 


A 

124 

A 

125 

A 

126. 

A 

127 

A 

128. 

A 

129 

A 

130 

A 

131 

A 

132 

A 

133 

A 

134 

A 

133' 

A 

13S 

A 

137 

A 

133 

A 

139 

A 

Uj 

A 

141 

A 

142 

A 

143 

A 

144 

A 

143 

A 

14&. 

A 

147 

A 

143 

A 

149 

A 

150 

A 

151 

A 

152 

A 

153 

A 

154 

A 

155' 

A 

156 

A 

157 

A 

153 

A 

159 

A 

160 

A 

161 

A 

162 

•  A< 

»163 

•A< 

»164 

A 

155' 

A 

165 

A 

167 

A 

163 

A 

169 

A 

170 

A 

171 

A 

172 

A 

173 

A 

174 

A 

173' 

A 

176 

A 

177 

A 

173 

A 

179 

A 

190 

A 

181 

A 

192 

A 

IB3 

A 

194 

A 

195 
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C  A  IBS. 

C        BALANCE  INJECTION  WELL  RATES  CONSIDERING  RESERVOIR  A  187 

C        THICKNESS  AND  STEADY  STATE  OPERATION  A  199 

C  A  199 

IF(NEDIT  ,GT.  0)  GO  TO  31  A  190 

SMQIsq.O  A  191 

SMQPso.O  A  192 

DO  20  I=1»NPIW  A  193 

IF  (Q(I).LT.O.O)  GO  TO  lQ  A  l9i 

SHQ1«SM0I*Q(I)  A  195' 

GO  TO  20  A  195 

10    CONTINUE  A  197 

SMQP»SMQP*Q(I)/HKI)  A  199 

20    CONTINUE  A  199 

DO  30  I=liNPIW  A  200 

IF  (Q(I).LT.O.O)  GO  TO  3o  A  201 

Q(l)a-Q(I)»SMQP*HI(I)/SMQI  A  202 

30  CONTINUE  A  20  3 
GO  TO  32  A  20V 

31  PRINT  595  A  205' 

32  CONTINUE  A  iOb 
C  A  207 
C  CHECK  IF  OIL  SANK  WILL  FORM  A  209 
C  A  209 

IF  (SOR.LT.SOI)  GO  TO  40  A  2lo 

PRINT  570  "  A  2n 

GO  TO  v70  A  212 

♦0    CONTINUE  A  213 

C  A  i\(> 

C        CONVERT  FIELD  UNITS  TO  COS  A  215 

C  A  21S 

DO  50  I*1»NPIBW  A  217 

XW(I)»XW(I)»30.*8  A  219 

YW(I)arW(I)ft30.*B  A  219 

QO(I)b9O(I)*1.8*01  A  220 

HI(I)bHI(I)»30.*8  A  221 

Q(I)eQ(I)*1.840l  A  222 

IF  (I.LEiNPIW)  30  TO  5o  A  223 

HIU)nH»30.*B  A  22* 

50    CONTINUE  A  225' 

DO  60  I*1»NBP  A  22«>. 

XB(I)aXB(I)»30.48  A  227 

YB(I)bYB(I)»30.*8  A  229 

60    CONTINUE  A  229 

H*H«3o.*8  A  230 

DO  70  I*1»NXH  A  231 

READ  590*  (HT(ItJ) tJsltNYH)  A  232 

70    CONTINUE  A  233 

READ  530.  (XH(I)»I=1.NXH)  A  234 

READ  580.  CYH(IJ  »I»1.NYH)  A  235' 

DO  80  I»ltNXH  A  236 

DO  80  U*1.NYH  A  237 

HT(I.U)=HT(I.J)*30.*8  A  239 

80        CONTINUE  A    239 

DO    90    IM.NXH  A    240 

XH(I)«XH(I)»30.4B  A    241 

90    CONTINUE  A  242 

DO  100  I»liNYH  A  243 

YH(I)bYH(I)»30.48  A  244 

100   CONTINUE  A  245' 

RI"RI»30.48  A  246 

DELT»DELT»B6400.  A  247 
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DELTT«DELTT*B6*00. 
DELPaDELP»86400« 
TMX=TMX*86400. 
0MNP«OMNP»liB*0l 

XMXbXMX»30.*B 

YMXsYMX«30;*8 

QMNl=QMNI»[.840l 
C 

C  SUBROUTINE  BOUND  CALCULATES  THE  NECESARY  RATES 

C        FOR  PREVIOUSLY  SELECTED  IMAGE  W^LLS  SO  THAT  THE 
C        STREAMLINES  ARE  CONFINED  WITHIN  THE  RESERVOIR 
C        BOUNDARIES. 
C 

CALL  BOUND 

PRINT  600 
C 

C  CHECK  QMNI  AND  CALCULATE  VELOCITIES  NEAR' 

C        INJECTION  AND  PRODUCTION  WELLS 
C 

DO  116  leltNPIB* 

IF  (Q(I).LT.O.O)  GO  TO  HO 

IF  (Q(I)/MI  (D.LT.3MNI)  QMNI sQ ( n /HI  ( I ) 
110   CONTINUE 

VLCIso.5»QMNI/RI 

VLCP*0»*#QMNP/H/RI 

NTSL»0 
C 
C 

C  THIS  SECTION  STARTS  A  STREAMLINE  ON  THE  CIRCLE  OF  A 

C        PRODUCTION  WELL  AND  FOLLOWS  IT  BACKWARDS  TO  LOCATE  THE 
C        PROPER  INJECTION  WELL.   ThIS  PROCESS  IS  CARRIED  OUT  FOR 
C        ALL  PRODUCTION  WELLS  AND  ALL  STREAMLINES.   THE 
C        STREAMLINE  STOPS  AT  A  POINT  WIThIN  THE  CIRCLE  OF  AN 
C        INJECTION  WELL.   THIS  POINT  BECOMES  THE  STARTING  POINT 
C        WHEN  THE  STREAMLINE  IS  FOLLOWED  IN  THE  FORWARD  OIPECTION 
C        TO  ACCUMULATE  THE>  PRODUCTION  HISTORY. 

C  AS  THE  STREAMLINES  ARE  BEING  FOLLOWED  A  PLOT  Opi  THEIR 

C        TRACE  IS  MADEt  (NOTE  THAT  ALL  PLOT  STATEMENTS  HAtfE 
C        ASTERISKS  In  THEIR  STATEMENT  NUMBERS.) 
C 

e 

C  PROGRAM  VARIABLES 

C 

c 

C  J  »  INDEX  FOR  WELLS 

C  K  «  INDEX  FOR  STREAMLINES 

C 

C  XIiYI  ■  POSITION  ON  STREAMLINE  OF  FLUID  PARTICLE  BEHNG  TRACED 

C  VXtVY  a  VELOCITY  OFi  FLUID  PARTICLE  IN  X  AND  Y  DIRECTIONS 

C  PP  *  PRESSURE  AT  PRODUCTION  WELL  (SINK) 

C  PI  *  PRESSURE  AT  INJECTION  WELL  (SOURCE) 

C  X(J»K) tY(J.K)  a  STARTING  POINTS  FqR  WATER  FRONT  FROM  INJ.  WELL 

C  XOUtK)  lYO(JtK)  «  STARTING  POINTS  FOR  OIL  FRONT  FROM  iNji  WELL 

C  TSUM(J.K)  a  SUMMATION  OF  TIME 

C  RGO(J.K)  a  STREAMLINE  GAS-OIL  RATIO 

C  OQ(JiK)  a  STREAMLINE  OIL  RATE 

C  CO(JtK)  b  STREAMLINE  CUMMULATIVE  OIL  PRODUCTION 

C  WQ(JtK)  b  STREAMLINE  WATER  RATE 

C  CW(J«K)  b  STREAMLINE  CUMMULATIVE  WATER  PRODJCTION 

C  SCL  «  PLOT  SCALE  FACTOR 

C  NSL  «  NUMBER  OF  STREAMLINES  FOR  THE  PRODUCTION  WELLI 


A 

2*9. 

A 

249 

A 

250 

A 

251 

A 

252 

A 

253 

A 

254 

A 

255' 

A 

25S 

A 

257 

A 

259 

A 

259 

A 

260 

A 

261 

A 

262 

A 

263 

A 

264 

A 

265' 

A 

266 

A 

267 

A 

269 

A 

269 

A 

27  8 

A 

271 

A 

272 

A 

273 

A 

274. 

A 

27  5' 

A 

275 

A 

27  7 

A 

279 

A 

279 

A 

230 

A 

281 

A 

282. 

A 

293 

•A«284 

•  A*283' 

•A< 

»296 

A 

28  7 

A 

289 

A 

289 

A 

290 

A 

291 

A 

292 

A 

293 

A 

294 

A 

295 

A 

29S. 

A 

297 

A 

299 

A 

299 

A 

300 

A 

301 

A 

302 

A 

303 

A 

304 

A 

305' 

A 

306 

A 

307 

•A 

»309 

A 

309 
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e  BEING  WORKED  ON 

C  OT  *  TIME  INCREMENT  FOR  SMALL  STEP  TAKEN.   STEP  SIZE  IS 

C  EQUAL  TO  .25**1.   RI  NEED  NOT  BE  SMALLER  THAN 

C  1/50TH  OF  THE  DISTANCE  BETWEEN  WELLS. 

C 

c 

c 

SCL*5.23B/YMX 

XM=2.25*XMX»SCL 

CALL  B3NPLT  (4LPLOT.XM.25.25) 

CALL  PLT  (1.125*XHX«sCLil.l25»YMX«SCLt-3> 

CALL  PLT  (XB(NBP)»SCLtYB(NBp)*SCL.3) 

DO  120  I=liN3P 

CALL  PLT  (XB(I)»SCL^YB(I)»SCLi2) 
120   CONTINUE 

DO  l*o  I»ltNPlw 

IF  (Q(I) .GE.0.0)  GO  TO  130 

CALL  SYMBOL"  <XW  ( I )  »SCL»  YW  ( I )  »SCLt  ,07»  1  lOfl ) 
GO  TO  1*0 
130   CONTINUE 

CALL  SYMBOL  (XW ( I ) »SCL» YW ( I ) »SCL« .07. 7»0»-l > 
1*0   CONTINUE 

DO  220  J*1*NPIW 

IF  (Q(J).GE.O.O)  GO  TO  220 

NSL=-NST«Q(J)/QMNP*o.5 
C 

C  THE  STREAMLINE  STARTING  POSITIONS  ARE  EVENLY  SPACED 

C        AROUMD  THE  PRODCTION  WELL 
C 

DO  210  Ksl.NSL 

XI«XW(J)»RI»COSU.0*K«6.283p/NSL) 

YI«YW(J)*RI»SIN(1.0*K«6.2B32/NSL) 

CALL  PLT  (XI»SCL»YI«SCLf3) 

PP(J.K)=0.0 

DO  150  L«lfNP!8W 

PP(J.K)=PP(J.K)-Q(Li)/HKL>#ALOG((XI-XW(L))»»2*(YI-YW(L))**2) 
150   CONTINUE 
160   CONTINUE 

IF  (ABS(XI).GT.XMX)  GO  To  180 
IF  (ABS(YI).GT.YMX)  GO  To  1B0 

VX»0.0 

VY=0,0 

DO  170  LaltNPIBW 

DEN=(XI-XW(L))*«2*(YI-YW(LM»*2 

VXsVX-0(L)»<XI-XW<L))/(DEN»Hl(L) ) 

VY«VY-3(L>#(YI-YW(L) ) / (DEN'Hl <L) ) 
170   CONTINUE 

VT«SQRT(VX»»2*VY**2J 

DT«0.25*RI/VT 

XI*XI«0T»VX 

YI«YIOT»VY 

CALL  PLT  (XI»SCL.YI»SCLt2) 

CHECK  IF  THE  STREAMLINE  HAS  REACHED  AN  INJECTION  WELL. 
THE  VELOCITY  OF  THE  FLUID  PARTICLE  IS  CHECKED  FIRST  TO 
SEE  I?  IT  IS  LAH3E  ENOUGH  TO  BE  IN  THE  VICINITY  OF  AN 
INJECTION  WELL.   IF  IT  IS,  EACH  INJECTIOM  WELL  IS  THEN  CHECKEO. 

IF  (VT.LE.VLCI)  GO  TO  l6o 
DO  200  LLBltNPlBW 

IF  (O(LL).LE.O.O)  GO  TO  200 
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310 
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311 

A 

312, 

A 

313 

A 

31* 

A 

315 

A 

31b 

A 

317 

•  A« 

313 

•  A* 

319 

•  A« 

320 

«A« 

321 

•  A« 

32? 

•  ft* 

323 

•A« 

32* 

«A« 

325' 

•  A« 

325 

*A« 

327 

•A« 

323 

»A« 

►3?? 

•  A« 

►330 

•  A* 

►331 

•  A* 

►332, 

A 

333 

A 

33* 

A 

335 

A 

335 

A 

337 

A 

339 

A 

33? 

A 

3*o 

A 

3*1 

A 

3*? 

•A»3*3 

A 

3** 

A 

3*5' 

A 

3*5 

A 

3*7 

A 

3*9 

A 

3*9 

A 

350 

A 

351 

A 

352 

A 

353 

A 

35* 

A 

355 

A 

35S 

A 

357 

A 

359 

A 

359 

A 

360 

A 

361 

•A 

»362 

A 

363 

A 

36* 

A 

365 

A 

366 

A 

367 

A 

369 

A 

369 

A 

370 

A 

371 
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RAD*SQRT({XI-XW(LL>>»»2»(Yl«.rirf<Lt)>**2>  *  372 

IF  (RAD.GT.RI)  SO  TO  200  A  373 

180   CONTINUE  A  374 

NTSLsNTSL*l  A  375 

X(JtK)eXI  A  374 

Y(JtK)«Yl  A  377 

TSUM(J,K)sO,0  A  373 

X0(JtK)«Xl  A  37? 

YO(JiK)«Yl  A  3S0 

PI(J«K)*0*0  A  381 

DO  190  LKsltNPIBW  A  382 

PI(JiK)=PI(J»K)-Q(LK)/HI(LK)»ALOG{  (XI-XW  (LO  )  »«2*  ( Y I-YW  UK)  )*«2»  A  38  3 

190   CONTINUE  A  394. 

RG0(U,K)=G0R(J)  A  385' 

IF  <SCW*S0I»SGR.GE.1,0>  RGO< J«K) *sGOR  A  386 

0Q(JtK)sQ0(J)/NSL  A  397 

CO(J.K)=o.O  A  399 

WQ(J«K>"0.0  A  39? 

GO  TO  210  A  390 

200   CONTINUE  A  391 

GO  TO  160  A  392 

210   CONTINUE  A  393 

CW(J.K)=0.Q  A  394 

220   CONTINUE  A  395' 

CALL  PIT  (0.0tOt0t?99)  »A»39S 

PRINT  610  A  397 

C  A  399 

C          END  OF  RUNBACK  SECTION  A  399 

C  A  400 

CTR"1.6  A  4fll 

PRINT  590  A  402 

TSTOPaDELT  A  403 

TSMPsOELT  A  404 

NBP1=N9P«1  A  405' 

230   CONTINUE  A  406 

SMCTYsO.O  A  407 

C  A  409 

C  A  4o? 

C          THIS  SECTION  STARTS  THE  STREAMLINE  AT  THE  INJECTION  A  4l3 

C        "ELL  POUND  ABOVE  AND  FOLLOWS  IT  IN  THE  FORWARD  DIRECTION  A  Ml 

C        TOWARD  THE  PRODUCTION  WELL.   THE  MOVEMENTS  OF  THE'  OIL  AND  A  41> 

C        WATER  FRONTS  ARE  CLOSELY  WATCHED  TO  DETERMINE  WHEN  THEY  A  413 

C        BREAKTHROUGH  INTO  THE  PRODUCTION  WELL.   THE  PRODUCTION  A  414 

C        HISTORIES  ARE  ACCUMULATED  AS  THE  FRONTS  MOVE.  A  415' 

C          LOOPS  420  AND  410  ARE  DONE  FOR  ALL  WELLS  AND  ALL  A  416. 

C        STREAMLINES  UNTIL'  TIME  EQUALS  TsTOP.   THEN  THE  CONSTANT  A  417 

C        PRESSURE  CORRECTION  FACTOR  <CTR)  IS  UPDATED  AND  THE  A  41 3. 

C        RESULTS  ARE  PRINTED  IF  DESIRED.  A  41? 

C  A  420 

C   »• ••»»#•••» ••• A  421 

C  A  4?2 

C               ADOITIONALi  PROGRAM  VARIABLES  A  423 

C  A  424. 

C      CTT  ■  CONDUCTIVITY  RATIO  A  425- 

C      CTR  *  CONDUCTIVITY  RATIO  CORRECTION  FACTOR  A  *2«, 

C      TSTOP  «  TIME  TO  STOP  AND  UPDATE  CTR  A  427 

t              TSMP  «  TIME" TO  PRIiT  AND  PLOT  A  429 

C      XIiYI  ■  LOCATION  OP'  WATER  FRONT  ON  STREAMLINE  BEING  FOLLOWED  A  42? 

C      XOIiYOI  «  LOCATION  OF  OIL  FRONT  ON  STREAMLINE  BEING  FOLLOWED  A  430 

C  A  431 

C   #••••••••• •••*••• • • ••* ••  A  432 

C  A  433 


00  42o  JsltNPIW 

IF  (Q(J).GE.O.O)  60  TO  420 

NSL*-NST»Q(J)/QMNP^o,5 

00  41o  KelfNSL 

CTY«0.0 

IF  <ABS(X(J,K)  )  .GT.XMX)  GO  TO  UO 

IF  (ABS(Y<U|K>  UGT.YMX)  GO  TO  4l0 

IF  (WQ(J»K).6T.0.OJ  GO  TO  400 

XI*X(J,K) 

YIsyjj.K) 

TSUMI»TSUM(J,K) 

XOl*XO(JtK) 

YOUYO(JtK) 

240 

c 
c 
c 

CONTINUE 

OETERMINE  THICKNESS  OF  RESERVOIR  A 

DO  250  MaltNXH 

IF  (XI.LT.XH(M))  GO  TO  260 

250 

CONTINUE 

MaNXH 

260 

CONTINUE 

DO  27o  N*1,NYH 

IF  (YI.LT.YH(N) )  GO  TO  280 

270 

CONTINUE 

NaNYH 

280 

CONTINUE 
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434 
435 
435 
437 
439 
439 
4«f) 
441 
442. 
443 
444 
445' 
445 
447 
449 
449 
450 
451 
45? 
453 
454 
455' 
455 
457 
459 
45? 
450 

HTP«HT(M,N)  A  461 

IF  (RGO(U,K) .LE.SGOR)  GO  TO  36o  A  462 

DO  29o  Msl.NXH  A  463 

IF  (XOI.LT.XM(M))  GO  TO  300  A  *54 

290   CONTINUE   "  A  465' 

MsNXH  A  466 

300   CONTINUE  A  467 

DO  310  N=1,NYH  A  469 

IF  (YOI.LT.YH(N))  GO  TO  320  A  469 

310   CONTINUE  A  47o 

N*NYM  A  4-71 

320   CONTINUE  A  472 

HTPO«HT(MiN)  A  473 

C  A  474 

C  THIS  SECTION  MOVES  THE  OIL  AND  WATER  FRONTS  AND  ACCUMULATES      A  475' 

C        PRODUCTION  BEFORE,  THE  OIL  FRONT  REACHES  A  PRODUCTION  WELL*         A  475 

C  A  477 

VXO«0.0  A  479 

VYO=0.0  A  479 

POB0.0  A  490 

VXaO.O  A  491 

VY«0.0  *  *B2 

PaO.O  A  493 

DO  336  L=1»NPIBW  A  494 

SQOs(XOl-XW(L) )##2*<Y0I-YW(L) )»»2  A  495 

P0«P0-3(L)/HI(L»»AL0G(SQ0)  A  495 

VX03VX0*Q<L)/HI (L) »(XOI-XW(L) »/SQO  A  497 

VY0*VY0*Q(L)/HKL)»(YOI-YW(L)  )/SQO  A  499 

SQ«(XI«Xm(L>  )»«2*<YI-YW(L)  )*»2  A  499 

P»P-Q(L)/Ml(L)»AL03(SO)  A  490 

VXaVX»3<L)/HI(UMXI-X*(U  ) /SO  A  491 

VYaVYOtU/HKUMYI-YrftL)  )/SQ  A  492 

330   CONTINUE  A  493 

IF  (PO.LT.PPU.K))  PO«PP(UtK)  A  494 

CTY«(PP(UiK)-PI(UtK))/(((PP(UiK)-PO)/8MB*(P0-P)/0M9*(P-PI(JtK) )/     A  495' 
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1       WMB)*GM8)  A  *95 

VLX*VX/(6.2832»P0R»(1-SCW»S0^-SGR) )»CTY/CT«»Hl(J)/HTP  A  *97 

VLY»VY/<6.2832»P0R*(l-SCW»SoR-SGR) > *CTY/CT«»Hl < J) /HTP  A  *99 

VELsSQRT(VLX»»2»VLY»»2)  A  *99 

DT»RI/VEL  A  500 

CRFB1.0*OQ<J.K)»NSL»CTR/(Q(J)»CTY)  A  501 

VLX0sVX0/(6.2832»P0R»(l-SCW-S0l-SGR)  ) »CTY/CTR»CRF»HI  (J)/HTPO  A  30? 

VLYO=vy0/(6;2832»P3R»(l-SC«'-SOI-SGR)  > »CTY/CTR°CRF»HI  < J)/HTPO  A  503 

VEL0»SQRT(VLX0»«2*VLY0»»2)  A  5q* 

IF  (VEL0.6T.VEL)  DT=RI/VEL0  A  503' 

IF  (TSJMIOT.GT.TSTOP)  DTwTSTOP-TsUMl  A  505 

XI*»XI»0T»VLX  A  5t)7 

YI«YIOT»VLY  A    509 

TSUMI=TSUMI*0T  A  509 

X0I=X0I*DT»VLX0  A  5lo 

Y0I»Y0I»0T»VLY0  A  5l 1 

C  A  512. 

C          CHECK  IF  OIL  FRONT  HAS  REACHEO  A  PRODUCTION  WELL.  A  513 

C        THE  WATER  FRONT  IS  ALSO  CHECKED  JUST  IN  CASE  IT  HAS  GOTTEN  A  5l  * 

C        AHEAD  OF'THE  OIL  FRONT.   THIS  OF  COURSE  IS  NOT  NORMAL  A  5l5> 

C        BUT  CAN  HAPPEN  I Fl  THE  OIL  PRODUCTION  RATE  DURING  A  515 

C        FILL-UP  IS  TO  LASGE.  A  517 

C  A  519 

IF  (SQRT((XI-XW(J) )»«2*(YI-YW(J))»*2).LE.RI)  GO  TO  3^0  A  519 

CO(JtK)=CO(JtK)*DT»OQ(JtK)  A  523 

IFJVELO  .LE.  VLCP)  GO  TO  35o  A  5?i 

IF(SQRT( (XOI-XW(J) )»»2MYOl-YW<j) )*«2).LE,RI)  GO  T0>  3*8  A  522 

DO  3*9  LX  =  liNPI8»(  A  523 

IF(U(LX)  .GE.  0.0)  GO  TO  3*9  A  52* 

IF(SQRT((XOI-XW(LX))»»2«(YOI-YW(LX))*#2).LE.RI)  GO  TO  348  A  525 

3*9   CONTINUE  A  52$ 

C  A  527 

C         WHEN  THE  OIL  FRONT  REACHES  A  pRODUCTlOM  WELL  THE  A  529. 

C       STREAMLINE  GAS-OIL  RATIO  IS  SET  EQUAL  TO  S60R.  A  529 

C  A  530 

GO  TO  350  A  531 

348   RGO(J,K)=SGOR  A  53? 

NCOUNT«l  A  533 

350   CONTINJE  A  53* 

IF  (TSUMI.LT.TSTOP)  GO  To  240  A  *35 

X(JtK)tiXI  A    535 

Y(J»K)bYI  A    537 

XO(J»K)=XOI  A    539 

YO(J»K)«YOI  A    53? 

TSUM<J.K)sTSjMI  A  5*0 

GO  TO  410  A  5*1 

C  A  5*2 

C         THIS  SECTION  MOVES  THE  WATER  FRONT  AND  ACCUMULATES  A  5*3 

C        THE  PRODUCTION  AFTER  OIL  FRONT  BREAKTHROUGH  A  5**. 

C  A  5*5 

360   CONTINJE  A  5*5 

VX"0.0  ft  5*7 

VYsO.O  *  5*9 

PeO.O  A  5*9 

DO  37o  L»1.NPI8W  A  5S0 

SQ»(XI-X»((L)  )**2«{YI-YW(L)  )»#2  A  551 

P*P-Q(L)/HI (L)*AL03«SQ)  A  552 

VX«VX*Q(L)/HI(L)*<XI-XW<L))/SQ  A  553 

VY«»VYO(L)/HI(L>»<ri-YW<L)  )/SQ  A  55* 

370   CONTINJE  A  555 

IF  (P.LT.PPUtK)  )  P«PP(JiK)  A  555' 

CTY«<PPU»K)-PI(Ji<n/<UPP<JtK)-p)/0M8*<P-PI(JtK>)/WMB)»GMB>  A  557 
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VLX"VX/(6.2832»POR»(l-SCW-So«-SGR) ) •CTY/CTRoHl U) /HTP  A  359 

VLVeVY/(6.2832*P0R»(l-SCrt-S0R-SGR))*CTY/CTR»Hl(J)/STP  A  55? 

VEL=S0RT(VLX»«»2*VLf»»2)  A  560 

DT»RI/VEL  A  5&i 

IF  (TSUMMDT.ST.TSTOP)  DT«TSTOP-TS<JMI  A  562 

XI«:XI»OT*VLX  A  563 

YI»YIOT»VlY  A  56* 

TSUMIbTSUMI*DT  A  565' 

OQ(J.K)=-Q(J)»CTY/(CTR»NSL)  A  56S 

CO{J»K)=CO(J,K)*OT*OQ<J.K)  A  567 

C  A  569 

C          CHECK  FOR  WATER  FRONT  BREAKTHROUGH  A  56? 

C  A  57Q 

IF  (SQRT((XI-XW(J))»«2*(YI-YW(J) )»»2).LE.RI)  SO  TO  390  A  37l 

3B0   CONTINUE  A  57? 

IF  (TSUMI.LT.TSrOP)  GO  To  2*0  A  573 

X(J»K)eXI  A  57* 

Y(JiK)=YI  A  373 

TSUM(U»K)=TSUMI  A  576 

XO(JiK)=XOI  A  577 

YO(JiK)=YOI  A  579 

GO  TO  *10  A  57? 

390   CONTINUE  A  58o 

C  A  581 

C          THIS  SECTION  IS  USED  WHEN  THE  WATER  FRONT  FIRST                  A  59? 

C        BREAKS  THROUGH,  A  593 

C  A  58V 

RGO(J, K)=SGOR  A  585 

CTY=WM9/6M8  A  586 

WQ(JtK)=-Q(J)»CTY/(NSL*CTR>  A  597 

OQ(J.K)=0.0  A  589 

X(JfK)aXI  A  59? 

Y(JtK)sYI  A  59Q 

TSUM(J,K)=TSUMI  A  591 

XO(J«.<)=XOI  A  592 

YO(J.K)=YOI  A  593 

CW(JtK)»(TSTOP-TSUMI)»WQ<J»K)  A  59* 

GO  TO  410  A  593 

♦  00   CONTINUE  A  596 

C  A  397 

C        ACCUMULATE  PRODUCTION  AFTER  WATER  BREAKTHROUGH  A  399 

C  A  59? 

CTY=WM9/GMB  A  600 

WO(J»K)=-Q(J)»CTY/(NSL#CTR)  *                                   A  601 

CW(J*K)=CW(UtK)»OELT*WQ(J(K)  A  602 

♦10   SMCTYaSMCTY*CTY  A  &03 

♦20   CONTINUE  A  &0* 

CTR«SMCTY/NTSL  A  605 

IF  (TSTOP.LT.TSHP)  GO  TO  *5o  A  60S 

C  A  607 

C          AT  SPECIFIED  TIME  INTERVALS  A  609 

C        ALL  STREAMLINES  ARE  SUMMED  AND  RESULTS  A«E  PRINTEO  A  bo? 

C        FOR  EACH  WELL'ANO.  FOR  THE  ENTIRE  FIELD.  A  Si  D 

C        PLOTS  ARE  MADE  OFi  THE  OIL  FRONT  AND  WATER  FRONT  POSITIONS.        »A»Sll 

C  A  612 

ORT»0.0  *  M3 

WRT=0.0  A  61* 

SOIL*0*0  A  615 

SWAT»o.O  A  MS 

00  ♦♦O  J»1»NPIW  A  617 

IF    (Q(J).GE.O.O)  GO  TO  440  A  619 

6R"0,0  A  61? 
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OR"0*o 
WR*0.0 
OPRD«0.0 
WPROsq.O 

NSL«-NST»Q(J)/QMNP*0.5 
DO  430  Kcl.NSL 
OPRD=OPR0*CO(JtK)/158990. 
WPR0»WPRD*CW(J.K)/158990. 
QR«6R*RG0( J»K)»OQ(J,K)/1.8*0l 
OR»OR*OQ(JfK)/l.B*01 
WRBHR*MQ(JtK)/l«6401 
♦30   CONTINUE 
GSR=6R/0R 
WOR=WR/OR 
EXpXW(J)/30.48 
WY»YM(J)/30t46 
PRINT  620»  JiEXtrtY 
PRINT  630t  GRtOK.WR: 
PRINT  6*0»  OPRDtrtPRDtGSRtWOR 
ORT=ORT*OR 
WRT=WRT*WR 
SOlLaSOlL^OPRD 
SWAT=SWAT*WPRO 
♦♦0   CONTINUE 

TMEI=TSTOP/86*oO.O 
PRINT  65o»  ORTtWRT 
PRINT  65lt  S0IL»SW4T 
PRINT  660*  TMEItCTR 
C 

C  MAXIMUM  PLOT  LENGTH  IS  U  FT. 

C        PLOT  LENGTH  IS  ESTIMATED  AMD  IF  IT  IS  GREATER  THAN 
C        11  FT»  THE  PLOTTING  IS  STOPPED. 
C 

IF(dl.O-TOTFT)  .LT.  PLTFT>  GO  TO  441 
TOTFT=TOTFT*PLTFT 
CALL  03PLT 
441   CONTINUE 

TSMPaTSMP*DELP 
450   CONTINUE 
C 

C  INCREMENT  THE  TO  NEXT  STOPPING  POINT. 

C        IF  TIME  IS  GREATER  THAN  TmX  -  STOP  CALCULATIONS. 
C        IF  OIL  FRONT  HAS  BROKEN  THROUGH  IN  AT  LEAST  ONE  STREAMLINE 
C        tNCOUNTeiiSET  DELT  *  DElTMO. 
C 

IF(NCOUNT.EQ.I)  DELT=OELTT 
TSTOP«TSTOP*D£LT 

IF  (TSTOP.GT.TMX)  GO  TO  460 
GO  TO  230 
460   CONTINUE 

CALL  EMDPLT 
470   CONTINUE 
480   FORMAT   (//t8Xt»NP  *  »tI3t*   NI  *  *tI3»»   *3»<  =  •lis**   ^Sp  s  #»I3 

It*   NST  »  »tI3> 
♦90   FORMAT   (//i8X.»H  «  •tFlo.St*   DELT  3  *»Fl0.5»#   THX  »  •♦Flo.l** 

1RI  =  •.F10.5) 
500   FORMAT   (/7t8Xt*DELP  «  »tFl0.5.«   XMX  *    •♦Flo.lt*   V^X  ■  *tFlo.l» 

1  •   QMMP  «  •tFlO.S) 
510   FORMAT   <//t8X.*POR  «  «>tF6.5*#   SC*  ■  »tF6.5t#   SOI  "  «tF6.5t*   SO 

1R  ■  •fF6.5t*   SGR  ■  »tF6.5) 
520   FORMAT   (//t 8X,»0IU  MOB  »  #tF10.5,#   "AT  M03'«  »tFlo.5t«   GAS  M03> 
U  *>tFl0.5> 


A 

620 

A 

621 

A 

622 

A 

623 

A 

624 

A 

625' 

A 

626, 
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630 
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631 
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632 
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633 

A 

63  4. 
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635' 
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636 
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637 
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639 
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63? 

A 

640 
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641 
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643 

A 

64  4 
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645' 

A 

646 

A 

647 

•A»643 

•A»64? 

•A»650 

*A«651 

»A*652 

•A*S53 

•A»654 

«A«553' 

•A*55S 
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657 

A 

653 
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65? 
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660 
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661 
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662 
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663 
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664 

A 

663' 
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665 
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667 
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66? 

•A*670 

A 

671 

A 
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674 

A 
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530 
540 
550 
560 
570 

560 
590 
595 
600 
610 
620 
630 

640 

650 
651 

660 
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20 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 

1  IS  NOT 
FORMAT 
FORMAT  ( 
FORMAT 
FORMAT 
FORMAT 
FORMAT  ( 
FORMAT 

1ATE  =  *t 
FORMAT  ( 

1  E10.3i* 
FORMAT  ( 
FORMAT  ( 

1E10.3) 
FORMAT  < 


(6F10.0) 

(6110) 

(6F10.2) 

(2F10.0) 

(•NO  OIL  BANK  WILL  FORM  BECAUSE  THE  INITIAL!  OIL  SATURATION 

GREATER  THAN  THE  RESIDUAL  OIL  SATURATION.*) 

(12F5.0) 

1H1) 

(8X.»INTER>JAL  WELL  BALANCE  DEFEATED*) 

(8X»»B0UNDING  COMPLETE*) 

(8X,»RUN3ACK  COMPLETE*) 

//»5X»*WELL,  NO  **I2»*  AT  »»E10.3.*  » *iElO. 3»*.*> 

(/♦10*i*GAS  RATE  s  *tElo.3»*.  OIL  RATE  «  *»F8.2t*.  WATER  R 

FB.2,».*) 

/»10Xi*CJM  OIL  *  *»E10.3»«.  CUM  WAT  *  *»£10.3»*.  GOR  ■  •• 

.  WOR  «  «»E10.3> 
/15X,*FIEL0,  OIL  RATE  *  #,El0.3,*t  WATER  RATE  »  *»El0.3) 
/15X,»FIELD  CUM  OIL  »  *,El0.3»*.  FIELD  CUM  WATER  *  #t 

/15X,*TIME  *  •f£lO«3#*DAYS«  CONDUCTIVITY  RATIO  «  *f£l0.3 


1  *.*t///) 

END 

SUBROUTINE  TO  BOUND  IRREGULAR  RESERVOIRS 
REF»  J.-K.  LIN,  PH.D.  OISSERTATlON*  THE  UNIVERSITY  OF 
TEXAS  AT  AUSTIN  (DECEMBER. 1972) • 

SUBROUTINE  SOUND 

COMMON/AAA/  X»  ( 1  25)  ,YW  ( 125)  .Hi  (125)  »(J(125)  .  XB  ( 125)  . 
1YB(125) ,NbP,NP.NI.X8W,H 

COMNON/BBU/  C (55,55) ,D(55) f 01 (55) 

DIMENSION  Xl(55)  ,ri  (55)  fXSB(lOO)  ,YSB(100)  .M3(100)  . 
1YIB ( loo) • AA( 100*55) ,BB( 100> «D5( 100) 

NBPlcN8P*l 

PAI»3, 14159 

L«NP*NI 

00  10  JMfNBtf 

Xl(J)sXM(J#L) 

YI(J)*YW fj#t> 

CONTINUE 

XB(NSP))«X8(1) 

Y8(N8Pl)sYB(l) 

CC»0.5 

DO  20  IslfNBP 

XH-0.5*(XB(I»l)*Xd(I) > 

VM*0.5*(YB(I»1)*Y8(JM 

DS(I>=S0RT((X6(I»1)-XB(1))»»2*(YB(I*1)-YB(I))**2> 

HDS*CC*DSU) 

TN"YB<I«1)-Y9(I) 

TD*XB(I«1)*X3(I) 

THETAsATAN2(TN»T0) 

BETA»THETA-0,5*PA1 

XIB ( I ) *XM«H0S*C0S (BETA) 

XS8 ( I ) eXM*  HOS*COS ( BETA  »P A I ) 

YIB(I)cYM*HOS*SlN(S£TA) 

YS8(I>«YM*HDS*S1N(BETA*PA1> 

CONTINUE 

mm* 

HsNdP 

DO  30  I»1*N8P 

DO  30  J*1»NBW 

FNUM»{XSB(I)-XI(J))«>*2*(YSB(I>-YI(J))**2 


A 

682 

A 

583 

A 

634 

A 

685 

A 

685 

A 

687 

A 

688 

A 

6B9 

A 

590 

A 

591 

A 

69? 

A 

593 

A 

59* 

A 

695' 

A 

59* 

A 

697 

A 

599 

A 

599 

A 

7oo 

A 

701 

A 

702 

A 

703 

B 

1 

8 

2 

8 

3 

8 

A 

8 

3 

8 

6 

8 

7 

8 

9 

B 

9 

8 

10 

8 

11 

8 

12 

B 

13 

8 

n 

8 

15 

B 

15 

8 

17 

8 

19 

8 

19 

8 

20 

8 

21 

8 

22 

8 

23 

8 

24 

8 

25 

8 

2b 

8 

27 

8 

29 

8 

29 

8 

30 

8 

31 

8 

32 

8 

33 

8 

34 

8 

3J< 

8 

36 

a 

37 

8 

39 

0 

39 

§ 

40 
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FDEN*(XlB(I>-XKJ)  )*>*2«(YI6(I)-YI(J)  >  *#2 
FACTaFNUM/^UEN 
AA(ItJ)»ALOG(FACT) 
30    CONTINUE 

DO  50  J»ltN8P 
SUM=0.0 
00  40  K=1.L 

FN«<Xl3<J)-Xrf(K) )*»2*(YIB(J)-YW(K) )  *>*2 
FOs(XS3(J)-XW(K) )**2*<YSB(J)-YW(K) ) *»2 
OSUM=Q(K)/HI (K)»ALOG(FN/FD) 
SUM=SUM*DSUM 
40    CONTINUE 

BB(J)=SUM 
50    CONTINUE 

DO  70  I31*N 
DO  70  K=ltl 
C(K»I)»0.0 
DO  60  J*1»M 

C<KiI>=C(KiI)*AA(J»I)«AA(J.K) 
60    CONTINUE 

C(l*K)xC(K.I) 
70    CONTINUE 

DO  »0  I=1»N 
0(1)20*0 
DO  60  J=liM 

D(I)=D(I)*dB(J)«AA(J,I) 
80    CONTINUE 

CALL  SOLVE(N) 
DO  90  J=l*Ntfw 
Q(J*L) =QI (J) »H 
90    CONTINUE 
RETURN 
END 
C 

C  SUBROUTINE  SOLVE  USED  By  SUBROUTINE  BOUND  »  SEE'  BOUND 

C        FOR  REFERENCE. 
C 

SUBROUTINE  SOLVE(N) 
COMMON/BBB/  C(55i55) iD(55) f 01 (55) 
DIMENSION  A(55«55) iIU(55) 
NN=N*1 
MsN 

DO  10  I  =  1»M 
A(I»NN)=D(I) 
ID ( I ) = I 
DO  10  J=1*N 
A(It J)sC(ItJ) 
10    CONTINUE 

K»l 
20    CONTINUE 
NR*K 
NC«K 

8  =  ABS(A  <K,iO  ) 
DO    40    I=KtM 
DO    40    J=K,N 

IF  (ABS(A(I.J) >-6)  40*40*30 
30    CONTINUE 
NR*I 
NC«J 

B"ABS(A(I»J)) 
40    CONTINUE 

IF  (NR-K)  70*70*50 


B 

41 

B 

*2 

B 

43 

B 

4*. 
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c 
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c 

15 

c 

15 
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17 

c 

19 

c 

1? 

c 

20 

»» 

21 

c 

?? 

c 

23 

c 

2* 

c 

23 

w 

25 

»» 

ta 

27 

c 

29 

29 
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50    CONTINUE  C   33 

DO  60  JaK»NN  C   31 

C"A(NR,J)  C   3? 

A(NR,J)=A(K,J)  C   33 

A(K,J)=C  C   3* 

60    CONTINUE  C   35' 

70    CONTINUE  C   35 

IF  (NC-K)  100,100.80  C   37 

80    CONTINUE  C   33 

DO  90  I=1»M  C   3? 

CaA(I,NC)  C   43 

A(I»NC)aA(ItK)  C   41 

A(I»K)=C  C   42 

90    CONTINUE  C   43 

IsIO(NC)  C   44 

ID(NC)=ID(K>  C   43 

IO(K)=I  C   45 

100   CONTINUE  C   47 

110   CONTINUE  C   43 

IF  (A(K,K))  120»180tl20  C   4? 

120   CONTINUE  C   53 

KK  =  tWl  C   51 

00  Uo  U=NK»NN  C   5? 

A(K»J)=A(K. J)/A(K»<)  C   53 

DO  140  I=1»M  C   54 

IF  (K-I)  130«1<*0,130  C   55 

130   CONTINUE  Z      55 

A(It J) sA(I,J)-A(It K)»A(K.J)  Z       57 

140   CONTINUE  C   53 

K=KK  Z      59 

IF  (K-N)  20tll0»150  Z      S3 

150   CONTINJE  C   51 

DO  17u  1  =  1. N  C   62 

DO  170  U=1»N  C   S3 

IF  (ID(J)-I)  173.160*170  C   54 

160   CONTINJE  C   55 

Uld)sA(JtMN)  C   55 

170   CONTINUE  C   57 

60  TO  1*0  C   53 

180   CONTINUE  Z      59 

PRINT  200  C   73 

190   CONTINUE  C   71 

RETURN  C   7? 

200   FORMAT   (19H  NO  UNlgUE  SOLUTION  )  C   73 

END  C   74 

C  *P*   1 

C       this  subroutine  plots  the  locations  of  the  oil  and  #p»  i 

C        WATER  FRONTS  AS  THEY  DEVElOPE.  ALL  DATA  IS  TAKEN  FROM            *>p»   3 

C        THE  MAIN  PROGRAM.  «p*   4 

C  #P*   5' 

SUBROUTINE  OBPLT  «P»   5 

COMMON/AAA/  XW(125>  .YWU25)  .Hi  (125>  »Q<125)  .X3U25)  .  «p*   7 

1YBU25)  fNdP.NP.NI.NBw.H  «P*   3 

COMMON/CCC/  X(80.20) .Y(80. 23) »XO(B0»20) .YO(90»20) •  #P»   9 

1WQ(B0.20) »SCL»XMX.YMX.NST.UMNP.RGO(80»20) *  *P*  10 

2S60R.NPIW.TMEI  *P*  11 

CALL  SYMBOL  ( 0.5, 1 3.5. . 1*. 7HTIME  =  .27o.0.7)                        *>P*  12 

CALL  NUMBER  (999.G »999.0» . H . TM£I ,270.  0  .0)  *p*  13 

CALL  SYMBOL  (999.0.999. 0 t . 14 t5H  DAYS. 270. 0 » 5)  #P*  U 

CALL  PLT  (1.125»XMX»SCL.l.l25»YMX»SCL»-3)  *P«  15 

CALL  PLT  (XB(NBP)»SCL»Yfa(NdP)»SCL,3)f  •?*    15 

DO  10  I3liN8P  *?•    17 
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CALL  PIT  (XB(I)«SCL.YB(I)»SCL*2)  •?•    13 

10    CONTINUE  *P*  19 

00  30  I=1»NPIW  *P*  20 

IF  (Q(I).GE.O.O)  GO  TO  2o  •?*    21 

CALL  SYMBOL" (XW(I)»SCL»YW(I)»SCLt. 07*1. 0»-l)  *P*  22 

60  TO  30  »P»  23 

20    CONTINUE  *P*  2* 

CALL  SYMBOL  (XW ( I ) *SCL» YW ( I ) »SCL» .07* 7.0»-l )  *P*  2, 

30    CONTINUE  *P«  2S 

DO  50  J=1»NPIW  »P*  27 

IF  (Q(J).SE.O.O)  GO  TO  5o  «P«  23 

NSLC-MST«0(J)/QMNP»0.5  *P*  2? 

DO  50  <=ltNSt  *P*  30 

IF  (A8S(X0(J»K) ) .GT.XMX)  30  TO  40  *P*  31 

IF  (ABS(YOUiK)  )  .GT.YMX)  GO  TO  40  *P*  32 

IF  (RGO(J,K) .LE.SGOR)  GO  TO  40  *P*  33 

CALL  SYMBOL  {XO <U» *) *SCL» YO (U*K) »sCL« . 07i31 ? 0»-l>  »P*  3* 

♦0    CONTINUE  «p»  33 

IF  (ABS(X(J,K)).GT.XMX)  GO  TO  50  *p»  3b 

IF  (ABS(Y(J.K) ) .GT.YMX)  GO  TO  50  •?*    37 

IF  (WQ(J,K) .GT.0.0)  GO  To  50  *P*  39 

CALL  SYMBOL  (X ( U* K) «SCL* Y (U. <) «SC|_» .07.40* Ot -1 >  »P*  39 

50    CONTINUE  «p«  *0 

CALL  PLT  (0.0»0»0t999)  *P«  *1 

RETURN  *P*  *? 

END  »P*  43 
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